 
%
clear all
close all


  disp('CLT DEMO')
  lw = 2; 
  set(0, 'DefaultAxesFontSize', 16);
  fs = 15;
  msize = 5;
  
rand('state',0)
%-------------------------------
figure(1)
sample = poissrnd(0.5, [1,300]);
hist(sample, [-1,0,1,2,3,4,5,6])
axis([-1 6 0 200])
print -deps 'C:\NPBook\ProbBasics\CLTfig1.eps'
%
figure(2)
rand('state',0)
S_300 = [ ];
for i = 1:5000
S_300 = [S_300 sum(poissrnd(0.5,[1,300]))];
end
hist(S_300, 30)
axis([100 200 0 500])
print -deps 'C:\NPBook\ProbBasics\CLTfig2.eps'